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ABSTRACT 

Context. Repeated observations of exoplanet transits allow us to refine the planetary parameters and probe them for any time depen- 
dent variations. In particular deviations of the period from a strictly linear ephemeris, transit timing variations (TTVs), can indicate 
the presence of additional bodies in the planetary system. 

Aims. Our goal was to reexamine the largely unstudied OGLE2-TR-L9 system with high cadence, multi-color photometry in order to 
refine the planetary parameters and probe the system for TTVs. 

Methods. We observed five full transits of OGLE2-TR-L9 with the GROND instrument at the ESO/MPG 2.2 m telescope at La Silla 
Observatory. GROND is a multichannel imager that allowed us to gather simultaneous light curves in the g', r', i', and z' filters. 
Results. From our analysis we find that the semi-major axis and the inclination differ from the previously published values. With the 
newly observed transits, we were able to refine the ephemeris to 2454492.80008(±0.00014) + 2.48553417(±6.4 x 10"^)£. The newly 
derived parameters are a = 0.0418 ± 0.0015 AU, r,, = 1.67 ± 0.05 Rj, and inc = 82.47° ± 0.12, differing significantly in a and inc 
from the previously published values. Within our data, we find indications for TTVs. 
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1. Introduction 

In the study of extrasolar planets, the observation of planetary 
transits has a prominent position. It allows us to determine sev- 
eral parameters that are not accessible by other means, shedding 
more light on the nature of the transiting planet. From the transit 
light curve itself, the planetary radius and the orbital inclination 
can be found by investigating its shape and depth. Together with 
radial velocity measurements, the knowledge of the inclination 
allows us to determine the true mass of the planet, which can 
then be used in conjunction with the radius to find the planet's 
mean density, allowing its composition to be constrained. With 
spectroscopic observations during transit and secondary eclipse, 
the chemical composition of planetary atmospheres can be stud- 
ied (e.g. Charbonneau et al. 2002 Sw ain et al.|2008 1 . 

Despite its advantages and increasing success, the transit 
method presents limitations for small and farther out planets: it 
is strongly biased towards close-in planets since the geometrical 
probability of a planet to transit decreases with orbital separa- 
tion and the transit signal becomes much less distinct for small 
planets, reaching down to a transit depth of only 0.01% for an 
Earth-sized planet orbiting a solar-type star. With increasing or- 
bital periods, scheduling becomes a major constraint on the dis- 
covery of planets, as transit signals become sparser in time, so 
long-term observations are required. 

The transit timing method has the potential to detect small 
transiting or non transiting objects such as moons, Trojans, or 
additional planets that can be difficult to find with current meth- 
ods. It makes use of the fact that the gravitational impact of these 
objects can be measured as dynamical interactions within the 
system, causing the time between successive transits to vary, and 
in turn the transit occurs earlier or later than expected ( |Holman| 



& Murray|2005 1. The transit timing method aims at discovering 
additional objects by searching for changes in the mid-transit 
times of close-in transiting exoplanets. Naturally, this technique 
is more sensitive to massive perturbers but it can be sensitive 
enough to detect planets of only a few Earth masses which are 
located near mean motion resonances with the transiting giant. 
Several searches for TTVs in exoplanetary systems have been 
undertaken up to now (e.g. Alonso et al.|2009 Diaz et al.|2008 



Holman et al. 2006 1. While no additional objects have been 



found with this technique, it has been used to pose limits on the 



existence of such objects in several cases (e.g. Rabus et al. 2009 



Adams et al. 2010 1. To detect TTVs one needs to observe a large 



number of transits and build up a database of mid-transit times 
covering several seasons. Additionally, the repeated observation 
of transits enables us to refine the planetary parameters. 

The subject of this paper, OGLE2-TR-L9 b, has been identi- 
fied as a planetary candidate in the publicly available data of the 
OGLE-II project (Udalski et al.|1997j by |Snellen et al.| ( |2007] i. 
It was recently confirmed as a transiting exoplanet with the ob- 
servation of a full transit and radial velocity measurements by 
Snellen et al. (2009 ). OGLE2-TR-L9 b h as amassof4.5±1.5 My 
and a period of approximately 2.5 days ( [Snellen et al.|2009] l or- 
biting an F3V star, at this point in time the hottest star known to 
host a transiting planet. 

In this paper, we present five new transits of OGLE2- 
TR-L9 observed with the GROND instrument ( |Greiner et al.| 
2008) mounted at the ESO/MPG 2.2 m telescope at La Silla 
Observatory during April and May 2009. We recorded each tran- 
sit simultaneously in four optical channels, g', r', i', and z'. 
Using this newly available data and reanalyzing the data ob- 
tained by |Snellen et al.|p009| l, we redefine the planetary param- 
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eters and conduct a first search for TTVs. An overview of the 
observations and the data reduction methods used will be given 
in Section l2] while a description of the data analysis and error 
estimation is presented in Section [3] Section |4] summarizes the 
obtained results. 



2. Observations and Data Reduction 

2.1. Observations 

We observed five transits of OGLE2-TR-L9 using the GROND 
(Gamma Ray Burst Optical and Near-Infrared Detector) instru- 
ment ( |Greiner et al.|200 8 ) mounted at the ESO/MPG 2.2 m tele- 
scope at La Silla Observatory. Designed with the aim of observ- 
ing gamma-ray-burst afterglows, GROND simultaneously ob- 
serves in 4 optical (g' r' i' z') and 3 infrared (J H K) channels 
allowing us to gather 4 optical light curves of each transit event. 
In the infrared channels, series of exposures of 10 sec each were 
obtained. Since dithering was switched off for a large part of 
the observations to improve time sampling and photometric pre- 
cision, no stacking of the infrared frames is possible with the 
standard way of background subtraction. We therefore ignore the 
JHK channels in our analysis. 

The observations took place during the nights of April 10, 
15, 20, 25 and May 15 2009, corresponding to the Epochs 177, 



179, 181, 183 and 191 based on the ephemeris given byjSnellen 



et al. ( 2009 1. During each night, we observed the full transit plus 
at least 20 min of baseline before ingress and after egress. For all 
observations, the exposure time was kept fixed at 46.4 s. During 
the night of April 10, we used the slow (~ 45 s) readout mode 
achieving a cadence of 2.05 min. We incurred some problems 
with guiding during the night of April 15, leading to unstable 
data quality. During this night, we were able to test the fast 
(~ 15 s) readout mode and observed with a cadence of 1.95 min. 
For the nights of April 20, 25 and May 15, we used the fast read- 
out mode, increasing the cadence to 1.31 min for April 20 and 
1.25 min for April 25 and May 15. 

In our analysis, we included the transit observed with 
GROND on January 27, 2008 by [Snellen et"aL] ( |2009l l. The ob- 
servation procedure was essentially the same as in our first ob- 
servation (April 10). A detailed description of the observations 
can be found in [Snellen et al.| ( [2009 ), while the data reduction 
was redone for this work as described below. 



2.2. Data Reduction 

With one of the goals of the initial OGLE observations being the 
detection of micro lensing events, OGLE2-TR-L9 is located in a 
rather crowded field. Fortunately, OGLE2-TR-L9 is reasonably 
bright (1=13.97 mag) and fairly isolated, with only very minor 
sources in its vicinity allowing the use of aperture photometry. 

The image calibration was performed with the mupipe soft- 
ware which has been developed at the University Observatory 
Municq^ and encompassed overscan and flatfield corrections. 
Then we carried out flux measurements for aperture photome- 
try using IRAPrldaophot procedures. We iteratively selected ap- 
proximately 15 of the most stable stars in the field (with the exact 



http://www.usm.uni- 



miipipe is available from 

[muenchen.de/people/arri/mupipe/ 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. The phased light-curves of OGLE2-TR-L9. From top to 
bottom, the filters are: g', r\ /', and z'. The phased light curves 
also contain the new reduction of the data from Snellen (2009). 
The best fit models are depicted as continuous lines. 



number depending on filter and seeing conditions) and combined 
them to a reference source to be used for differential photometry. 
For comparison, we also performed aperture photometry us- 
ing mupipe but found that the light curves produced with IRAF 
show a smaller scatter Therefore we only used the light curves 
produced with IRAF for the analysis. Comparing the light curves 
in the four optical channels, we found that the r' filter shows the 
best accuracy, with a photometric accuracy of 1.5 mmag, while 
we achieved 1.8 mmag for the g' channel, and 2.0 mmag for the 
i' and z' channels. We present all light curves in Figurel6] 

3. Determination of System Parameters 

The light curve fits were performed using the model for a plan- 
etary transit given by Mandel & Agol (2002). We assumed a 
quadratic limb darkening law and chose the coefficients accord- 
ing to Claret (2004) for a star with metallicity [FejH] - 0.0, sur- 
face gravity log{g) - 4.5 and effective temperature T^ff - 7000 
K. We found the best fit to our data by minimizing the x' using 
a downhill simplex algorithm as implemented in the AMOEBA 
code ( [Press et al.|[1992| l. Our goal was to find the best-fitting 
values for the central transit time t^, the ratio of the planetary 
and stellar radii r;,/r,, the planetary semi-major axis in units of 
the stellar radius alr^, and the orbital inclination inc. The above 
parameters have the advantage that they can be derived directly 
from the light curve without assumptions regarding stellar prop- 
erties other than limb darkening. As starting values, we assumed 
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Fig. 2. The results for the radius ratio rp/r* obtained from the 
phased hght curves against the filter wavelengths. The values 
show a slight trend towards smaller planetary radii with decreas- 
ing wavelength. 



the central transit times calculated according to the ephemeris 
given by Snellen et al., (2009 ), and typical values of planetary 



transits for the other parameters: rplr^ = 0.1, a/r* - 20 and 
inc = 90°. After the best parameters were found for each light 
curve, we used the resulting mid-transit times to phase fold the 
data producing the light curves shown in Figure [T] To improve 
on the parameters derived from the data in each of the four fil- 
ters, these phase folded light curves were again fitted for r^/r,, 
a/r^ and inc. The final values were found by a combined fit of 
all available light curves. The resulting parameters are shown 
together with the results for the phased light curves in Table [T] 
and the corresponding model is depicted as a continuous line in 
Figure [T] 

To estimate the errors of the determined parameters, we used 
the Bootstrap Monte-Carlo Method (Press et al. 1992) which 
works by creating a large number of representations of the data 
by randomly choosing, with replacement, subsets from the orig- 
inal data set. This means that each newly created data set con- 
tains the same number of points as the original data set but with 
some points left out and some points duplicated. For each boot- 
strap data set, the best fitting parameters are found using a pro- 
cedure identical to the analysis of original data and the errors 
are calculated from the distribution of the results. To get a sec- 
ond estimate for the errors in the central transit time, we kept 
the model parameters fixed at their best values and let f^ vary 
until the deviation in the x^ exceeded Ax^ - 1 . The resulting 
errors are in good agreement with the errors derived from the 
Bootstrap Monte-Carlo Method and range from 29 seconds for 
transits observed with a lower cadence to 16 seconds for transits 
observed with higher cadence. Here, it should be pointed out that 
the above error treatment does not incorporate the contribution 
of correlated low frequency noise ( Pont et aL]|2006) . Thus, the 
errors on the mid-transit points are likely to be underestimated. 



4. Results 

4.1. Planetary Parameters 

In Table [T] we present the results for the phased g', r', i', and 
z' light curves. As can be seen in Figure l2] the values for the 
planetary radius mainly agree within one sigma, apart from the 
value derived from the g' filter which lies slightly below. While 
values for the planetary semi-major axis and inclination found 
with the g', r' and i' photometry agree within one sigma, the 
values derived from the z' band observations favor a larger or- 
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Fig. 3. The results for the planetary semi-major axis (upper 
panel) and the inclination (lower panel) against the filter wave- 
lengths. The fact that the shapes of the variations are very similar 
shows the degeneracy present between the semi-major axis and 
the inclination: larger solutions for the planetary separation favor 
higher inclinations. 



Filter: 


g' 


r' 


i' 


z' 


all 


rp/r. 


0.11021 


0.1116 


0.1121 


0.1123 


0.1120 




±0.0009 


±0.0006 


±0.0006 


±0.0006 


0.003 


air. 


5.99 


5.80 


5.93 


6.23 


5.88 




±0.15 


±0.11 


±0.11 


±0.13 


±0.06 


inc [°] 


82.69 


82.31 


82.56 


83.25 


82.47 




±0.32 


±0.25 


±0.24 


±0.26 


±0.12 



Table 1. The results for the planetary to stellar radius ratio, the 
semi-major axis in units of the stellar radius, and the inclination 
for the filters used next to the final parameters produced by a 
combined fit of all available light curves. 



bital separation combined with a higher inclination (see Figure 

As pointed out by Seager & Mallen-Ornelas (20031), the 
mean stellar density can be derived directly from the light curve 



shape and the planetary 



Seager & Mallen-Ornelas 



period. Following the description of 
( 2003) , we find the stellar density to 
, which is in good agreement with the 



bep, =0.442 + 0.014 [po 
previously published value. 

To find the planetary parameters Up and r^, we adopted the 
stellar parameters from Snellen et al. (2009) and combined them 
with the best-fit values obtained from our analysis. The results 
and their respective errors can be found in Table l2| All values 
are in good agreement with the values published in S nellen et al.| 
(12009) except for Up and inc where our new results are signif- 
icantly different. The reason for this discrepancy is an error in 



the calculation of fl„ in Snellen et al. (2009i where 1/M» was 
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Parameter 



Semi-major axis [AU] 

Planetary radius [Rjup] 

Inclination [°] 

Transit duration [d] 

Period [d] 

To [HJD] 

Stellar density [pg] 



0.0418 ±0.0015 
1.67 ±0.05 
82.47 ±0.12 
0.1092 ±0.0021 
2.48553417 ± 6.4 x 10"' 
2454492.80008 ± 0.00014 
0.442 ±0.014 



Table 2. The final results for the planetary parameters derived 
from the data presented in this work. 



t, [HJD - 2400000] 


Reference 






50478.661 ± 0.0012 


Snellen etal. 2007| 


54492.80086 ± 0.00033 


ttiis worJi 


54932.73944 ±0.00031 


this work 


54937.70976 ± 0.00029 


this work 


54942.68169 ±0.00022 


this work 


54947.65232 ±0.00019 


this work 


54967.53689 ± 0.00022 


this work 



Table 3. All mid-transit times of OGLE2-TR-L9 known up to 
date. 



used instead of M* when transforming the measured period into 
a semi-major axis using Kepler's third Law. The error in semi- 
major axis propagated into a wrong inclination value. We recom- 
mend to use the values pubUshed here for future work because 
of this calculation eiTor in .Snellen et al.| ( |2009j and because the 
new values are based on a significantly larger dataset. 



4.2. A New Ephemeris 

To derive a new ephemeris for OGLE2-TR-L9 b, we used all six 
observed transits of OGLE2-TR-L9 together with an ephemeris 
derived from the OGLE-II dataset by [Snellen etaT (2007). It 
should be noted that this "mid-transit point" is derived from an 
OGLE-II light curve which was constructed from ~ 500 data 
points with irregular cadence collected over a period of 3.5 years. 
From our new analysis of the transit observed by [Snellen et al. 
( 2009| l on January 27, 2008, we find the central transit time to be 
HJD 2454492.80086 ±0.00033 instead of HJD 2454492.79765 + 
0.00039. All mid-transit times known for OGLE2-TR-L9 can be 
found in Table |3] Fitting a constant period to all available mid- 
transit times gives the following new ephemeris; 



Tc(E)[HJD] = 2454492.80008(+0.00014) 
H-2.48553417(±6.4 x 10"^)£. 



-2000 -1500 -1000 -500 
Epoch 




Epoch 




Epoch 



Fig. 4. The O-C diagram for the ephemeris calculated from all 
known mid-transit times of OGLE2-TR-L9 b. In the upper panel, 
all points are shown, while the lower two panels zoom in on the 
new points together with the point found from the transit in 2008 
(middle panel) and solely the new points (lower panel). The mid- 
transit points of the epochs and 179 occur 1.5 min later and 1.0 
min earlier than expected, respectively. 



correspond to the newly observed transits. It can be seen that the 
transit times show some deviations from a constant period as the 
transit of epoch occurs 1 .5 min later and the transit of epoch 
179 occurs 1.0 min earlier than e xpec ted. Fitting a constant pe- 
riod to the data as done in Section 
a confidence level of > 99.75%. 



4.2 



gives ^;^^ = 3.86 yielding 
However, one must take into 
account that the errors on the central transit times are likely to 
be underestimated due to correlated noise (see Section[3]l. 

Next to variations in the transit timing, we also searched our 
data for variations in the transit duration which could hint at the 
existence of exo-moons (Kipping 2009p. For this purpose, we 
measured the transit durations based on the best models of both 
the phased light curves and the light curves derived from each 
transit event. The results are illustrated in Figure Island show no 
detectable variation. 



4.3. Transit Timing and Transit Duration Variations 

Finally, it is possible to search for TTVs. To do so, the observed 
minus calculated (O-C) values with reference to the ephemeris 
determined above are computed for each of the given mid-transit 
points. The result is shown in Figure [4| Here epoch -1615 cor- 
responds to the ephemeris derived by Snellen et al. (2007) from 
the OGLE-II data, epoch corresponds to the corrected mid- 
transit time from |Snellen et al.[ ( |2009| l and the last five epochs 



5. Conclusion 

We have gathered five full transits of OGLE2-TR-L9 b with the 
GROND multichannel imager mounted on the ESO/MPG 2.2 m 
telescope at La Silla Observatory. With this new data, and in- 
cluding the transit observed by S nellen et al.| ( 2009] l in our anal- 
ysis, we recalculated the parameters of OGLE2-TR-L9 b and 
found different semi-major axis and inclination values compared 
to [Snellen et al.[ ( |2009| l. This is due to a calculation error in the 
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Fig. 5. The transit durations obtained from the phased Hght 
curves in the four fihers (upper panel) and the Hght curves for 
each transit event (lower panel). All values agree within their 
errors. 

previous work and thus the values presented in this work should 
be used in the future. We studied the central transit times and 
transit durations for any variations that could be attributed to a 
perturbing body in the OGLE2-TR-L9 system. While the tran- 
sit durations agree for all transits, the mid-transit points show 
indications of period variations. 
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Fig. 6. All light curves in the four optical channels. The data were recorded during the nights of January 28, 2008 (a) (by^ 

let al.|2009| l, April 10 (b), April 15 (c), April 20 (d), April 25 (e) 2009 and May 15 (f) 2009. The continuous lines represent the best 
fit to the respective data. 



